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ABSTRACT 



We examine recent developments in the cluster cooling flow scenario following 
recent observations by Chandra and XMM-Newton. We show that the distribu- 
tion of gas emissivity verses temperature determined by XMM-Newton gratings 
observations demonstrates that the central gas, i.e., where the cooling time is 
less than the age of the cluster, in cooling flow clusters cannot be in simple 
steady-state, i.e., M is not a constant at all temperatures. Based on the mea- 
sured gas emissivity, the gas can only be in steady-state if there exists a steady 
heating mechanism that scales as H{T) oc T° where a = 1 — 2. That is, a heat- 
ing mechanism that preferentially targets the hottest and highest entropy gas, 
which seems very unlikely. Combining this result with the lack of spectroscopic 
evidence for gas below one-third of the ambient cluster temperature is strong ev- 
idence that the gas is heated intermittently. While the old steady-state isobaric 
cooling flow model is incompatible with recent observations, a "moderate cool- 
ing flow model" , in which the gas undergoes intermittent heating that effectively 
reduces the age of a cooling flow is consistent with observations. Most of the gas 
within cooling flows resides in the hottest gas, which is prevented from cooling 
continuously and attaining a steady-state configuration. This results in a mass 
cooling rate that decreases with decreasing temperature, with a much lower mass 
cooling rate at the lowest temperatures. Such a temperature dependent M is re- 
quired by the XMM-Newton RGS data and will produce an increasing amount of 
intermediate temperature gas which will then be reheated during the next heat- 
ing cycle. We show the compatibility of this model for the cooling flow cluster 
A2052. The present paper strengthens the moderate cooling flow model, which 
can accommodate the unique activities observed in cooling flow clusters. 

Subject headings: galaxies: clusters: general — cooling flows — intergalactic 
medium — X-rays: galaxies: clusters 
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1. Introduction 

Several recent papers have shown that the predictions of the steady-state cluster cooling 
flow (CF) model are inconsistent with Chandra and XMM-Newton X-ray observations (see 
review by Fabian 2003). Chandra observations show that the gas in the central regions 
of relaxed clusters with central dominant galaxies can be described as homogeneous single 
temperature gas with a positive temperature gradient. Ettori (2002) showed that the ASCA 
evidence for multiphase CFs was due to the poor spatial resolution of ASCA which could 
not distinguish between multiphase gas and single phase gas with a temperatiu^e gradient. 
Chandra observations show that only within the central few tens of kpc in clusters does 
the spectroscopy require additional components above a single temperature. Of course, this 
could simply be due to statistical and spatial resolution limitations in the Chandra data. 
The strongest evidence against the steady-state isobaric CF model is the lack of observered 
line emission in XMM-Newton RGS spectra from gas cooler than ~ 1/3 of ambient cluster 
temperatures (e.g., Kaastra et al. 2001; Peterson et al. 2001; Fabian 2003). In the so called 
"standard steady-state CF model" , the age of a CF is assumed to be similar to the age of 
the cluster; namely the gas has been cooling for a long time. In the so called "moderate CF 
model" (Soker et al. 2001) heating is intermittent and the hot gas is not in steady-state, in 
the sense that its cooling time is longer than the time elapsed since the gas was last heated. 
Many reheating scenarios of the central gas in CFs have been proposed in the past with 
the aim of suppressing CFs altogether, or significantly reducing the average mass cooling 
rate (e.g., Binney & Tabor 1995; Tucker & David 1997; Ciotti & Ostriker 2001; David et al. 
2001; Quihs, Bower, & Balogh 2001; Briiggen, M. & Kaiser 2001; Ruszkowski & Begelman 
2002; Nulsen et al. 2002). These scenarios can be divided into models with steady heating 
by conduction or AGN, or non-steady, self-regulating models heated by nuclear outbursts. 
The moderate CF model is in the class of non-steady models where the gas cooling rate 
is peridocially heated by nuclear outbursts. We present below the evidence for non-steady 
CFs. The frequent occurrence of X-ray cavities coincident with radio lobes around the central 
dominant galaxy in CFs (e.g., Abell 496, Dupke & White 2002; Perseus [Abell 426], Fabian et 
al. 2002; Hydra A, McNamara et al. 2000) demonstrates that AGNs have a significant impact 
on the X-ray morphology of the hot gas. However, there are still significant uncertainties in 
the detailed physics of how the relativistic and thermal plasmas interact and how much heat 
is deposited into the hot gas. 

In the moderate CF model where the intermittent heating is generated by AGN activity 
in the central dominant galaxy, it is expected that the amount of cooling gas increases sharply 
with increasing temperature. Only at very low temperatures is there a steady-state situation 
(i.e., where the cooling time is short compared to the time between nuclear outbursts). In 
Soker et al. (2001) we presented the arguments for a moderate CF in which the actual mass 
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cooling rate is significantly below that derived under the assumption of the standard model, 
but a non-steady CF still exists. Soker et al. (2001) estimate that in the moderate cluster 
CF model the required kinetic energy of the AGN is ~ 10^^ erg s~^, and its strong activity 
should last ~ 10'^ yr and occur every ~ 10^ yr. Only ~ 1% of all CF clusters should be found 
during that stage. In Cygnus A there is a strong radio source which heats the intra cluster 
medium (Smith et al. 2002). Wilson, Young, & Smith (2003) estimate the mechanical power 
of the jets in Cygnus A to be L^et — 6 x 10^^ erg s"-*^. This is ~ 100 times larger than the 
radio emission, and it is in the range required by the moderate CF model. 

In recent analyses of XMM-Ncwton RGS data by Kahn et al. (2003) and Peterson et al. 
(2003; hereafter P2003) they find that there is less gas than that predicted by the steady- 
state isobaric CF model at all temperatures below the ambient gas temperature. Also, the 
discrepancy increases with decreasing temperature. There is not just a deficit in gas below 
~ 1/3 of the ambient cluster temperature, there is a deficit of gas at all temperatures below 
the ambient temperature. They conclude that their results are difficult to reconcile with 
the newly proposed alternatives to the standard CF model, including those that completely 
suppress radiative cooling with some form of steady-state heating, and that new physics may 
be required. Our goal in this paper is to show that the XMM-Newton results are consistent 
with the expectations of the moderate CF model presented in Soker et al.(2001). 

We convert the differential luminosity as a function of temperature as derived from the 
RGS data on cluster CFs into the distribution of gas mass verses temperature in §2. We 
also derive M as a function of temperature in this section and show that the gas in CFs 
cannot be in steady-sate. In section 3 we apply the moderate CF model to a recent Chandra 
observation of A2052, and in §4 we summarize our main results. 



2. Distribution of Gas Mass with Temperature 

We show in this and the next sections that the distribution of gas mass with temperature 
within cluster CFs found by P2003 can be incorporated into a model with intermittent 
heating, such that the effective age of the cooling gas is only ~ 1 — 3 x 10^ yr. 

Based on RGS spectra of 14 CF clusters, P2003 find that the variation in the differential 
luminosity with gas temperature can be characterized by the following expression: 

dT 2 iinip \To J 

where k,/!, and nip have their usual meaning, Mocf is the inferred mass cooling rate based on 
the assumptions inherent in the steady-state isobaric CF model, T is the gas temperature. 
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and To is the maximum, or ambient cluster temperature. In general, the luminosity of gas 
cooling isobarically within a temperature interval dT is: 

In steady-state, the rate gas cools into a given temperature interval must equal the rate 
gas cools out of the same temperature interval. In other words, M{T) must be a constant 
and equal to Mocf- Comparing equations (1) and (2) shows that the gas can only be in 
steady-state if a = 0. Fitting the RGS data on their sample of 14 CF clusters, P2003 find 
that q; ~ 1 — 2. Within the old CF radius— where cooling time equals the cluster age— most 
of the radiation comes in the X-ray band. Only in the very inner region of r < 10 — 30 kpc 
a significant fraction of the energy lost by the cooling gas may be emitted in the optical and 
UV band (e.g., Fabian et al. 2002; Soker, Blanton & Sarazin 2003). 

The mass cooling rate as a function of temperature can be written as: 

M(r, ^ l^f ^ 1^-^, (3, 

al at al Tcooi 



where we define the cooling time to be 

dT/dt 



Tcool(r) = (4) 



^From equations (2) and (3) we obtain 

dL 5 k dM T 



dT 2 fxrup dT TcooI 
Combining equations (1) and (5) gives: 



(5) 



dM ■ , /T\"~' 1 

MoefTcool(«+l)h^ 7^. (6) 



dT ' \ToJ Tc 

The cooling time varies as: 

AoTno AoPo/tV 

where tq is the cooling time of gas at Tq, n is the total number density, P is the gas pressure, 
and A is the radiative cooling function, such that An^ gives the energy radiated per unit 
volume per unit time. The last two equations can be combined to give: 

MocMa + 1) — 777^-^- (8) 



dT ' V ^0/ To A P 



- 5 - 



The cooling function can be characterized as: 

A /r\^ 



A \T I 

i-O \ / 

where ^ 1/2 for T > 2 x 10^ and ~ — 1/2 at lower temperatures. For consistency with 
the expressions derived above, assuming isobaric cooling and integrating equation (8) over 
temperature gives: 

a + 1 / T \ """^^ 
M(< T) = MocfTo ^ - (10) 



a -77 + 2 \To^ 

Setting a ~ 1.5 in the last equation, the average value found by P2003, and assuming 
Bremsstrahlung cooling (r) — 0.5), gives: 



rj-i \ 3 



M(< T) ~ MocfTo (^-J . (11) 

If we add a heating mechanism to equation (2), then a steady-state condition can be 
estabhshed only if H{T) oc T". It is difficult to conceive of a heating mechanism that prefer- 
entially heats the hottest and highest entropy gas in a CF. The RGS data show that M{T) 
increases with increasing temperature. If this were true over the lifetime of a cluster, there 
would be a large reservoir of gas at intermediate temperatures, which is also inconsistent 
with the RGS data. This gas must be periodically removed from these intermediate tem- 
peratures either by cooling sporadically to very low temperatures, which simply returns us 
to the classic CF problem, i.e., the lack of a significant reservoir of cool gas, or intermittent 
heating back to roughly the ambient temperature. We therefore examine the non-steady 
moderate CF model. 



3. Moderate Cooling Flow Model 

As an example, we consider the CF cluster A2052, which was included in the P2003 
XMM-Newton sample, and whose X-ray structure as observed by Chandra was discussed in 
detail by Blanton et al. (2001, 2003). P2003 assume a cooling radius of ro = 51 arcsec, and 
find q; ~ 3 for this cluster. Hence, from equation (10), M(< T) oc T^-^, and most of the 
mass is in the hottest gas. 

Indeed, from fig. 4 of Blanton et al. (2003) we find that the mass within 30 arcsec 
consists of ~ 30% of the total gas mass within 50 arcsec. The temperature of the gas at 30 
arcsec is 2.5 keV (Blanton et al. 2003). Using this along with an ambient temperature of 
kTo = 3.3 keV (Blanton et al. 2003), imphes that (2.5/3.3)^-5 = 29% of the total gas mass 
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up to To resides at temperatures lower than 2.5 keV. Considering the uncertainties, these 
two numbers are in excellent agreement. A small fraction of the gas still resides at lower 
temperatures, presumably because the shock that heated the gas, say ~ 10^ yrs ago, could 
not increase the cooling time of the lowest entropy gas above the time between outbursts 
(Soker et al. 2001). 

P2003 did not compare their XMM-Newton results directly with Chandra data. We 
find the cooling time at 30 arcsec, with kT = 2.5 keV and rig = 0.02 cm^"^ (Blanton et al. 
2001), to be rcooi(30kpc) = 1.5 x 10^ yr. Interior to this radius the intracluster medium is 
disturbed by two large radio bubbles (Blanton ct al. 2001). In the moderate CF model, 
the time interval between intermittent energy deposition is 1 — 4 x 10^ yr, hence the gas is 
prohibited from settling into a steady-state configuration at r > 30 arcsec. Although the 
gas continues to cool and its cooling time gets shorter, it is unable to reach a steady-state 
before the next nuclear outburst. 

Only within r ~ 30 arcsec may the gas have reached a steady-state, however the Chandra 
data show that it was recently disrupted by the two radio bubbles. We therefore expect that 
the actual coohng rate is much smaller than that in the old (standard) CF model. Based on 
a spectral analysis of the Chandra data, Blanton et al. (2003) find the mass coohng rate to 

be 26 < M < A2Mq yr-\ which is ~ 1/3 of the old value within r ~ 140 arcsec (Pcrres et 
al. 1998). Taking the actual cooling radius to be the radius within which steady-state has 
been established based on the shorter age in the moderate CF model, the cooling rate will 
be even lower than the value found by Blanton et al (2003), i.e., we argue for a mass cooling 
rate of M < lOM© yr^^ 

Although we only study one cluster as an example, we note: (1) In the moderate CF 
model suggested by Soker et al. (2001), the heating doesn't inhibit the CF in the very inner 
region r < 10 kpc. In this region, the gas continues to cool to temperatures of T ~ 10"^ K, 
even after a heating event. Only in the outer regions does the heating event prevent the gas 
from cooling to low temperatures. Therefore, we don't expect the gas to be isothermal. The 
heating event can't heat the gas to extremely high temperature either, because this requires 
AGN energy output much larger than typical observed values. (2) The moderate CF model 
thus predicts that a in eq. (1), or a — + 2 in equation (10) will not take extreme values. 
We can't predict the exact range of values of a in the present study; this requires numerical 
simulations with variable conditions in the intracluster medium before each heating event, 
the energy supplied by the event, the time elapsed between events, and other processes, e.g., 
mergers. (3) Some CF clusters have central cooling times shorter than those in A2052. This 
does not pose a problem for the moderate CF model since the central regions of clusters 
still harbor a CF. (4) In light of the uncertainties, e.g., the temperature profile in the inner 
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regions, the presence of X-ray cavities, and in the model parameters mentioned in point (2) 
above, it does not warrant a more extensive comparison with other clusters at this point. 
Future work, in particular numerical simulations of heating events, will include a greater 
comparison between the moderate CF theory and cluster observations. 

4. Summary 

We show that the differential luminosity of the gas in cluster coohng flows as a function of 
temperature, as derived from RGS XMM-Ncwton observations, is inconsistent with steady- 
state cooling flow scenarios, but is consistent with non-steady heating or moderate cooling 
flow models. The findings of P2003 imply that most of the gas in cooling flows resides in 
the highest temperature phase (cqs. 10 and 11 above). Within the context of the moderate 
cooling flow scenario (Soker et al. 2001), at these temperatures and densities the CF cannot 
reach a steady-state, and if the intermittent heating continues, it cannot attain such a 
state. Therefore, the rate of gas coohng at high temperatures is much higher than at lower 
temperatures. 

We demonstrate the applicability of this model to the cooling flow cluster A2052 (Sec. 
3). We argue for a mass cooling rate of < IOMq yr~^ in this cluster, which is compatible 
with the upper hmit found by P2003 for lowest temperature gas. This is much lower than the 
coohng rate of hotter gas, > lOOM© yr~^, found by (P2003) or deduced in the old CF model 
by Peres et al. (1998), and somewhat lower than the cooling rate derived recently by Blanton 
et al. (2003) of ~ 26 — 42Mq yr~^. Overall, intermittent heating in the moderate CF model 
(Soker 2001; Fabian 2003), or more frequent heating ( Blanton et al. 2003 for A2052), may 
account for the P2003 findings without invoking new processes. The intermittent heating 
model, with time intervals between major heating events of ~ 1 — 3 x 10^ yr has the advantage 
that no fine tuning is required to balance heating and cooling, since the gas is heated to 
relatively high temperatures, and then starts coohng. Most of the heated gas does not cool 
to low temperatures before the next major heating event. Different values of the physical 
parameters, e.g., energy input and time intervals between heating events, as well as a cluster's 
properties, will give different values of a in equations (1) and (10). Indeed, a large range of 
values is observed, a ~ 1 — 3, hinting on a wide range in the physical parameters mentioned 
above, such that no fine tuning is observed or required. As shown in Soker et al. (2001), it is 
difficult to prevent the lowest entropy gas from complete cooling (i.e., increasing the cooling 
time above the time between heating events). Hence, a low M cooling fiow can be sustained 
in the very central region of clusters. 

This paper strengthens the moderate cooling fiow model (Soker et al. 2001), by support- 
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ing the claim that the unique activities observed in coohng flow clusters (e.g., McNamara 
2002 and references therein), can be accommodated within its framework. The problems 
which are need to be solved, e.g., the exact nature of the heating events, are less severe than 
the crisis encountered in the old cooling flow model. 

We thank an anonymous referee for useful comments. This research was supported in 
part by grants from the US-Israel Binational Science Foundation and G02-3171. 
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